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Abstract: We analyze the possibility of primordial magnetic field amplification by 
a stochastic large scale kinematic dynamo during reheating. We consider a charged 
scalar field minimally coupled to gravity. During inflation this field is assumed to be 
in its vacuum state. At the transition to reheating the state of the field changes to a 
many particle/anti-particle state. We characterize that state as a fluid flow of zero 
mean velocity but with a stochastic velocity field. We compute the scale-dependent 
Reynolds number Re(k), and the characteristic times for decay of turbulence, td and 
pair annihilation t a , finding t a <C We calculate the rms value of the kinetic helicity 
of the flow over a scale C and show that it does not vanish. We use this result to 
estimate the amplification factor of a seed field from the stochastic kinematic dynamo 
equations. Although this effect is weak, it shows that the evolution of the cosmic 
magnetic field from reheating to galaxy formation may well be more complex than 
as dictated by simple flux freezing. 
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1. Introduction 



The question of the origin of large scale magnetic fields that permeate almost all 
structures of the universe is one of the most challenging areas of research in astro- 
physics. None of the main lines of investigation, namely primordial origin or in situ 
generation, succeeded up to now to explain both the intensity and the topology of 
the large scale fields. Local generation mechanisms are mainly based on seed field 
generation by, e.g., a local battery, amplified by a turbulent dynamo in the interstel- 
lar or intergalactic medium (see |IJ and references therein). The primordial origin 
hypothesis, on the other hand, considers that at least the seed field is generated 
at some early epoch (inflation, reheating or radiation dominance), and is amplified 
by flux conservation and/or turbulent dynamo action during gravitational collapse 
from z ~ 100 on [Q. The seed field must be quite intense for gravitational collapse 
to produce the detected intensities, and the turbulent dynamo must operate almost 
since the birth of the galaxy, i.e., during most of the matter dominated era. The 
recent detection of regular fields in high redshift quasars 0, @, || however may 
challenge the in situ generation, or at least the dynamo mechanism in the form we 
understand it today, favoring the primordial origin of the fields. 

Two obstacles must be overcome by a successful primordial generation mecha- 
nism: breaking the conformal symmetry of a massless gauge field in a spatially flat 
universe and building a large coherence length. Sub-horizon processes, like phase 
transitions 0, 0, 0, 0, in general produce intense fields, but of very small co- 
herence length (see Refs QTOJ and |TI| for reviews of different magnetogenesis mech- 
anisms). The inflationary epoch of the universe (if ever existed) offers a suitable 
scenario for large scale field generation as in it sub-horizon scales naturally become 
super-horizon. Several mechanisms were considered along the years in which confor- 
mal invariance is broken either by coupling the magnetic field to curvature in different 
ways or addressing non-linear electrodynamics [0, fl~3 |, [[14], ||15|| . In general the 



fields produced are extremely weak, or of marginal intensity, to seed subsequent am- 
plification processes. The reheating period has also been studied as a magnetogenesis 



scenario 



17|, JIJ, pi, 



but in all scenarios considered so far the obtained 
fields are too weak to be of astrophysical interest. 

Confronted with this situation one wonders if it is possible to have a pre- 
amplification (or perhaps full amplification) of a seed field created by one of the 
above mentioned mechanisms already in the early universe. In this sense the reheat- 
ing epoch offers a good prospect, as it is a period where highly non-linear and out of 
equilibrium processes take place |21|, |22], |3|, |4], ^| 



91 30, Bll B2L This 



possibility was explored for the first time some years ago by Finelli and Gruppuso 
[3~3|j an d by Bassett et al. ||34j| . In Ref. |33] it is analyzed the amplification of a 



pre-existing magnetic field by parametric resonance during the oscillatory regime of 
a scalar field to which the magnetic field is coupled. In Ref. MM the amplification 
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during preheating is studied considering several different models. Another possibility 
for such pre-amplification process, and that will be investigated in this paper, could 
be the operation of a turbulent large scale dynamo P3| , p6f , JT[, |ft7| , similar to the 
one that acts in the interstellar plasma. 

That the matter fields in reheating can be turbulent was pointed out in Refs. 



[ PD| ] (see |HJ for a theoretical analysis of turbulent reheating). A 
dynamo requires the presence of a plasma. As the inflaton is a gauge singlet, it will 
not decay directly into charged species. Therefore to have a plasma we must consider 
an extra, charged, field. The mechanism by which the plasma is created is particle 



creation during the transition from inflation to reheating |39|, |2T], [40 1. 

Suppose that the charged species in question was in its vacuum state during 
inflation. The created particles will generate stochastic currents that on one hand 
induce a seed field [ 16|| , [|18| and on the other may constitute the turbulent plasma 
we are looking for. Creation of spin 1/2 particles such as electrons is suppressed 
by conformal invariance at the high energies prevailing during inflation [IB], so the 
charged species must be a scalar. Suitable candidates can be found in supersymmetric 



extensions of the standard model 17 



The simplest model for a turbulent large scale dynamo is driven by flow velocities 
and does not take into account the back-reaction of the amplified fields. It is known 
as a kinematic dynamo |35|, [ 36[] , [0- The sufficient condition for it to be 



operational is the flow to be helical, i.e., that the volume average of the scalar product 
of the vorticity (curl of the velocity) and the velocity, known as kinetic helicity [[11|] , 
does not vanish [Q, [|43[]. Of course this approximation (the neglect of the back 
reaction of the induced field) is valid for weak magnetic fields and/or very short times 
of operation. Mathematically speaking, the equation for that dynamo can be written 
as dB l /dt ~ —t corr 'H c EijkdB k /dxi ) where B l is the large scale field (or mean field), 
% c the kinetic helicity and t corr a correlation time. VLe^dB^ jdx^ ~ B l /L, with L the 
coherence length of the field, then we can estimate B l (t) ~ Bq (0) exp {—t corr H c ti 1 L). 

In this paper we shall investigate the possibility of a dynamo action during 
reheating. We assume the existence of a charged scalar field, minimally coupled to 
gravity, that is in its vacuum state during inflation. To simplify the analysis we 
consider de Sitter inflation and thus a de Sitter invariant vacuum for the field |0] . 
As mentioned above, when the transition from inflation to reheating takes place, the 
scalar field is amplified, and stochastic currents are generated. The characterization 
of these particles as a fluid is straightforward. The hydrodynamic energy and pressure 
are determined by matching the expectation value of the energy-momentum tensor 
of the scalar field to that of a perfect fluid at rest. 

The fluid has a stochastic Gaussian velocity, which is found by matching the self- 
correlation of the Oi components of the energy momentum tensor of the fluid to the 
symmetric expectation value of the corresponding operator for the field. Finally, the 
viscosity of the fluid is found by assuming that it is close to saturate the Kovtun, Son 
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and Starinets bound [45]. While initially derived from consideration of the AdS/CFT 
correspondence, the fact that a similar bound seems to hold for the strongly coupled 
quark gluon plasma ]46| suggests that this bound is a good description of field theories 
in general. 

We characterize the turbulence by finding the momentum dependent Reynolds 
number Re (k). As for the magnetic Reynolds number, R m we do not need to estimate 
it because we are interested in the kinematic regime, where magnetic fields are too 
weak to backreact on the flow. As there are no stirring forces, turbulence will decay 
eventually. We calculate the decay time of the turbulence for each mode, t d (h). 
On the other hand, the fluid is made of particles and antiparticles, which are liable 
to annihilate. We also estimate the characteristic time for pair annihilation, t a (k), 
finding that t a {k) < t d (k), i.e., the fluid annihilates before turbulence decays. This 
fact allows us to consider that the turbulence is stationary in the interval < t < t a . 

The non-trivial result of our paper is that the rms value of the kinetic helicity of 
the fluid, He, is not zero. This proves that the kinematic dynamo action mentioned 
above is indeed possible. The key ingredient to have "H c 7^ 0, is the fact that the 
plasma is made up of two scalar fields, $ and & . We then estimate the amplification 
factor of the induced field based on the dynamo equation written above, i.e., ~ 
exp [H c t corr /L]. 

We work with signature (— , +, +, +) and with natural units, i.e., h = c = ks = 
1. Greek indexes denote space-time coordinates while latin indexes refer to spatial 
coordinates. To simplify our analysis we shall consider de Sitter inflation, and define 
dimensionless variables and fields as x l = Hr l , r = Ht, and $ = H~ l ty, where H is 
the Hubble constant that characterizes the de Sitter phase. The mass of the scalar 
field will combine with H to produce the dimensionless mass parameter, m/H. In 
section I we make a brief review of dynamo theory. Section II is devoted to the 
fluid description of a quantum field. In it we find the velocity correlation function 
and velocity spectrum as well as the kinetic helicity correlation function. In section 
III we characterize the turbulence by finding the Reynolds number, Re (k), and the 
characteristic times t d and t a . In section IV we find the amplification factor for the 
magnetic field. Finally in section V we summarize our conclusions. The bulk of the 
calculations that lead to these results are shown in the four appendices. 



2. Basics of mean field dynamo theory 

In this section we briefly sketch the so called first order smoothing approximation 
(FOSA) approach to the theory of mean field dynamo. We refer the reader to Refs 
[j35|| , HJ, ||7) and references therein for details. In FOSA, purely hydrodynamic 
turbulence is considered, ignoring higher than second order correlations in the fluc- 
tuating velocity field u l . This approach is suitable for short times and for magnetic 
fields that are weak enough to neglect their backreaction on the turbulent flow. In 
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short, fields are divided in mean and fluctuating components, as 



B l = B i 



U l = U i + v l 



(2-1) 



where overbar denotes volume average, and that is assumed that satisfies Reynolds 
rules [|^]. In the case that U l = 0, the mean magnetic field satisfies the equation 



dB i 



jjk oc 



- + 



dt dxi a 2 (r) 

The important quantity here is the mean electromotive force, S, given by 



£ l = tijkV j b k 

If B l is sufficiently weak and regular, £ l can be expanded as |l[ 



o L 



QB p (t') 

{t,t')BP(t') + (3 ikp {t,t') 



dx k 



dt' 



(2.2) 



(2.3) 



(2.4) 



with a ip = ei jk Vj (t) dv k (#) /dx p and /3 ikp (t,t') = e ikp vi (t) v p (f) . Under the hy- 
pothesis of local homogeneity and isotropy, the tensors in the integrand must be 
proportional to S ip and e ikp respectively eq. ( [2.4|) can be written as 



£' 



a {t- t') B { (t') +P(t- t') J* (t') 



dt' 



(2.5) 



where a (t — t') = — (1/3) v (t) ■ w (£'), w* = e 13 djU}~ being the vorticity of the velocity 
fluctuations; (3 (t — t') = (1/3) v (t) ■ v (f) and J 1 (t) the mean electric current. If 
besides it is assumed that B l (t) is a slowly varying function of time, then eq. ( |2.5| ) 
turns into 

£i = aB 1 - pji (2.6) 



with 



and 



a 



1 f l l 

-- / v(t) • Vl(t')dt « --WrV(t) -W(t) 

•J Jo <j 



/3 



1 



V (t) • v (t')dt' « -t corr v 2 (t) 



(2.7) 



(2.J 



with t corr . the correlation time. The last approximation in eqs. ( f2.7| ) and ( |2.8|) is 
known as the "r approximation" [|l]]. Observe that a is minus the kinetic helicity, 
T-L c of the flow |^T[. A non-null value of this quantity indicates that the flow lacks 
mirror symmetry. This is a sufficient condition for dynamo action |5B|, [37]. If B p 
is smoothly varying, then the dominant term in eq. (|2.6| ) is the first one, and eq. 
( |2. 2D can be written as 

dW 

(2-9) 



dt 



1 8B k 

3 ox 3 
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Taking €ijkdB k jdx^ ~ B l /L, with L the scale of coherence of the mean field, eq. 
Q2.9|) can be directly integrated for short times. We shall show below that in our 
case % c is a Gaussian variable of zero mean value and known variance £% . Taking 
the ensemble average over all possible realizations of H c we obtain that the mean 
magnetic field is 

W{t) -Sfexp Q(^Lr%^ 2 )) (2-10) 

with B l the initial value of the field. Our task in the next sections is to characterize 
the system of cosmological created scalar particles as a turbulent flow, and investigate 
if it has a non-zero kinetic helicity. 



3. Fluid description of charged quantum scalar fields 



Consider a charged scalar field, , minimally coupled to gravity in a spatially 

flat Friedmann-Robertson- Walker universe, described by the line element dS 2 = 
—dt 2 + a 2 (t) (dx 2 + dy 2 + dz 2 ), with a (t) the expansion factor. We assume that 
the e.m. field is so weak that it can be neglected throughout. The action of the field 
reads 



S 



-— I d 4 Xy 

2 



-g 



g a/3 d a $dp& + ^$$t 
H 



(3.1; 



with gv v the spacetime metric, m/H the dimensionless mass parameter of the field 
and H = a (r) /a (r) the Hubble constant during inflation. Throughout the paper 
we consider m/H -C 1 (see e.g. [0). The stress energy tensor is given by 



Explicitly 



H 4 



rpflU 



-2 55 



-9' 



-g>»'d && x & - -g^^r^ 
2 2 i7 2 



The electric current density is 



(3.2) 



(3.3) 



(3.4) 



We identify with the stress energy tensor of a two fluid system. One fluid 
corresponds to the positively charged scalar particles, and the other to the negatively 
charged anti-particles. Analogously we identify J£ with the electric current of the 
two fluid system. To this purpose we define the four velocity of the flow as usual, 



i.e., 



=j(U" + v") (3.5) 

with = (1,0,0,0) the four velocity of the fiducial observers at rest with respect 
to the radiation field produced by the decay of the inflaton. We define the projector 
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onto the surface orthogonal to the world lines of fiducial observers in the usual way, 
i.e., = + UW, so h^U v = and hTv v = v". We then write 

n o = (p + p) 1 2 U U°+pg 00 (3.6) 

Tf } = (p + p) 7 2 f/V, vWV + v*. (3.7) 

4 = en 7 (vV - vL) (3.8) 

where we have symmetrized Tj ^ = (T^ 1 + T^) /2. In eq. ( |3.7| ) fvM is the 
stochastic velocity of the positively (negatively) charged species and in eq. (|3.8|) n 
is the number density of particles. In both equations 7 is the Lorentz factor due 
to the (macroscopic) velocity of the fluid measured by the fiducial observers. Our 
flow is made up of gravitationally created particles during the transition between 
inflation and reheating. As momentum is conserved in the particle creation process, 
in the radiation frame both fluids have zero bulk velocity, thus we can take 7 = 1. 
Therefore v!j_ are stochastic fluctuations around the zero mean velocity, that must 
be characterized through their correlation function. 

3.1 Transition from inflation to reheating: particle creation 

The stochastic velocity v l introduced at the beginning of this section is the result 
of random motions of scalar charges. To understand how those charges appear we 
observe that a state detected as vacuum by inflationary observers will be detected as 
a many-particle state by comoving observers in the reheating epoch. Mathematically 



this is expressed as follows [21, 38, 39, 40 



From eq. (|3.1|) we obtain the evolution equation for the scalar field $, the Klein- 
Gordon equation, 

d 2 2 m 2 



dT t - V- - W ) * = (3-9) 

(and an identical equation for $'). The field operators can be written in terms of 
creation and annihilation operators as 

1 C Hk 

$ = 77^ J -T(7y [** W ^ + « W ( 3 - 10 ) 
$t = 7^72- / -T777 [*« ^ e ^ + « ^ e ~^4] (3-11) 



(2*)' 



where k = k/H, is the dimensionless wavenumber, k the physical wavenumber and H 
the Hubble constant during inflation. Replacing in eq. (|3.9|) we obtain the evolution 
equation for each mode <f> K , i.e., 



K 2 ( m \ 2 ^ (® -|_ ^ 



a 2 (r) \HJ 2\a 2a 2 



0. (3.12) 
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Here k/cl{t) is the dimensionless physical wavenumber. For simplicity we consider 
de Sitter inflation, where an invariant vacuum for a minimally coupled scalar field 
exists 0] , and we assume that the scalar field is initially in this state. Therefore the 
positive energy solutions of eq. ( p. 12 ) are the Hankel functions 



MB8I 



HP 



K 



a(r) 



(3.13) 



with v = ^9/A-m 2 /H 2 . We follow Refs. f4|, Q to the effect that during the 
reheating period the scale factor of the Universe scales as t 2//3 . In this case it is 
accurate enough to consider a WKB approximation for the modes, i.e., 

exp [-iS K (r)] 



tR 



T 



(3.14) 



with dS K (r) /dr = Vt K (r) = y K 2 /a 2 (r) + (m/H) . After the transition to reheating 
the commoving observers in the new geometry see the state of the scalar field as a 

Mathematically this is described as 

(3.15) 



many-particle state, [[38|, [39|, |2T|, |40[ 

<P{ (r) = a K <f>« (r) + /3 K 0f (r) 



where (j) T K (r) (0f (r)) are the positive frequency solution of the Klein Gordon equation 



for inflation (reheating), and a K and (3 K the so called Bogoliubov coefficients |38|, ^9], 
|. If /9 re 7^ 0, eq. ( p.!5|) shows that a positive frequency wave during inflation 



21 



becomes a mixture of positive and negative frequency waves during reheating. The 
details of the calculation of these coefficients are given in Appendix A, here we quote 
the resulting expressions together with the physical explanation. The number of 
created particles in modes with k < 1, i.e., super-horizon ones, is not sensitive to 
the details of the transition. For k > 1 that number does depend on the transition 
features. We take into account this dependence by assuming the most simple form 
for it, i.e., a linear transition that lasts a time tq. The (3 K coefficient for a linear 
transition reads 



(») 



9 

16 



1 1 1 



3 i2r 5[0] 



sin (2t q k) 



(3.16) 



where r is the duration of the transition from inflation to reheating, (see Ref. 
pT[] for a similar analysis, though for cosmological perturbations). As k > 1, the 
modes that are most amplified are those for which tqk is small. For these modes 
sin (2tqk) ~ 2roK. Hence we take 

9 \ 2 1 1 e i2ToS M 
16/ Ato k 5 

As for k < 1 details of the transition do not matter, we have the usual solution from 
assuming an instantaneous transition at r = 0. 

i(u-l)V(u-l) 1 



(3.17) 



(0 



7rV2fiV2 (Q) 



K 



(3.18) 
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with v = a/9/4 — m 2 /H 2 , and Q K (0) = a/k 2 + m 2 /H 2 . After the transition an out 
of equilibrium plasma is established. It has no bulk velocity with respect to the 
comoving observer's rest frame, but due to the random motions of its constituents, 
fluctuating velocities do exist. 



3.2 Two Point Velocity Correlation Function 

One way to characterize a system with fluctuating velocities is to give their spatial 
two point velocity correlation function fl52f . It is defined as the equal time ensemble 
average of the product v* (r, x) v J (r, x'), i.e., 



R 3 (r, x, x') = v* (r, x) v^' (r, x') 



(3.19) 



where supra-indexes indicate the Cartesian components of the turbulent velocity. 
From eq. (|3.7p, we can define a stochastic velocity operator as 



T. 



{Oi} 



T,X 



(P + P) (t) 



(3.20) 



and we assume that it is not relativistic. Observe that this does not mean that the 
particles are non-relativistic, they indeed are at such high energy. However their 
collective motion can be safely taken as non-relativistic. A state of a quantum field 
is specified by its Hadamard two point function, i.e., the vacuum expectation value 
of the anticommutator of the field at different spacetime points. So from definition 
( ggOp we can calculate ^0 |tJ 0i} (r, x) , T^ 0j} (r', x')} 
spatial two point correlation function of the velocity fie 



Oy, and from it obtain the 
d as 



R v (r, x, x) = lim 







Tf } (r,x),TrV,x 



{Oi} 







(p+py(r) 

using eqs. ( |3.10| ) and ( |3.11| ), the Hadamard two point function reads 



(3.21) 







Tf } (r,x),Tf' } (r',*')} 







II 



32vr 3 a 6 (r) J J 



dndw 



dr' 



, K i„j ® r i+ t- _n ^ r i+ f _/\ , .. ;^, i)2(: l' (j, T ') r iA 

^G ro (r, r j — G K (r, r J + G ro 

4 7 d 2 Gl+(r,T') 

+ K K - 



(3.22) 



dr'dr 



< + (r, r') 



with £ 



.r 



x', and where G{ H 



i>£ (r) 0^.* (r') is the positive frequency 



propagator. Writing the scalar field modes during inflation, 0^ (r) in terms of the 
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modes during reheating, 0^ (r) as (fij, (r) = (r) + (3 K (p^* (r) with a K and /3 K the 
Bogoliubov coefficients, we find the positive frequency propagator 



GI + (t,t') 



G« + (r, r') + (r) ^ (r') + (r) 0f , 

+ |/3 K | 2 [G^(r,r') + Gf-(r,rO] (3.23) 

with (r, r') = <fi^* (r) 0^ (r') the negative frequency propagator. When replac- 
ing eq. fl3.23| ) in eq. ( |3.22| ) there appear three kernels, one with vacuum contributions 
only, another with contributions both from the vacuum and from the created parti- 
cles, and a third one, built from contributions from the created particles alone. This 
is the most important one. Details of the calculations, as well as explanations of the 
approximations made along the way, are given in Appendix C. Here we quote the 
main results and discuss the physics involved. Replacing the propagators and their 
derivatives, and neglecting rapidly decaying terms, we obtain 

R ij (t,x,x + £)~ lim 00 _^ 8 fi/ _, ; 1 , 2 / / dhdw e l{R+ ^ \(3 K \ 2 |/3 ro | 2 (3.24) 



r'^r 327T 3 a 6 (t) (p + p)' 



w<w* [G* + (r, r') + Gl~ (r, r' 



+ kV[G? + (t,t')+G*- (t,t' 



d 2 



dr'dr 
d 2 



G« + (t,t>) + 



d 2 



G« + (t,t>) + 



dr'dr 
d 2 

-G 



G R K ~ (r, r') 



R- 



pr'dr ro v ' ' dr'dr 

The quantity (p + p) is calculated in Appendix B, and to obtain it we neglected its 
fluctuations. This means that we are identifying it with its expectation value. The 
result is 



(P + P) 



H 4 



r 2 + 



1 

24 



(3.25) 



2(2vr) 1/2 a 4 (r) 

Observe that it depends on two parameters, m/H and r which are related to the 
contribution of super-horizon and sub-horizon modes respectively. The velocity spec- 
trum, $ y (q, r), is given by the Fourier transform of eq. ( |3.24j ), i.e., 

1 



d 3 ^'(r,f,f + e)e- 



(3.26) 



As shown by Tomita et al 



(2tt) 3/2 

eddies larger than the horizon are frozen in the plasma 



and decay with the expansion. We shall consider only modes inside the particle 
horizon, i.e., modes that are in causal connection, so q > 1. Further calculations are 
sketched in Appendix C. The main contribution to the velocity spectrum is due to 
sub-horizon modes, almost parallel to q. After performing the calculations eq. fl3.26| ) 

reads 

-2 



& j r) 



x 



9 

16 



'TIT 



3 h 4 r 



512 vr H 3 



3 (l£ 1 
2 \m 



r 2 o + 



9 

16 



1 

24 



(3.27) 
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The general form of $* J (?, r) for isotropic turbulence can be written as [52] 



$ y {q, t) = [$ LL (s, r) - $ NN (s, t)} ^- + $ NN (s, r) S 13 



(3.28) 



where §ll {s, t ) is the longitudinal part of the spectrum and &nn 0?> t ) the normal 
part. By direct comparison of eq. (|3.28| ) with ( |3.27| ) we have 



®nn(s,t) ~ — — - 



9 \ 4 3 la 4 (r) / # x 4 



16 7 512 tt H 3 Vm ' T ° 



3 (H 



2 to' r ' 



9 

16/24 



- (3.29) 



9 y 3 l a 4 (r) /if 
16/ 128 7r if 3 V m 



'0 



3 (H 



2 V m 



9 

16/24 



(3.30) 



Observe that both functions are of similar amplitude. The energy spectrum is defined 

as 



9 V 9 a 2 (r) (H 



(3.31) 



16/ 128 H \m 



a /_9_ 
2 \m) T ' V 16 7 2 1 



-2 



where dfl (?) is the solid angle element, and 

$ 11 {s,t) = 2$ nn (?,r) + $ LL for) 



(3.32) 



The total energy per mass unit is then 



-] —a(r)(- 
16 J 256 w V ™ 



2 V m) T ° + V 16 J 24 



(3.33) 



We may define the total energy associated to a scale k > 1 as given by the contribu- 
tion of all eddies smaller than k~ 1 , i.e., 



E K 



1 



a it) 



■£(<T,t) 



(3.34) 



and so 



3.3 Kinetic helicity two point correlation function 



A sufficient condition to sustain large scale dynamo action is that the turbulence be 
helical |35|, |36], j37j . As stated in Sect. II, it is defined as the volume average of the 

i.e., 



scalar product of the vorticity and the velocity |41 



n. c . 



1 



WV 



Vol 



d (vol) wV 



(3.36) 



Vol 



with w J = e l i k djVk the vorticity of the velocity field. A non-null value of this quantity 
indicates that the flow lacks of mirror symmetry. Due to conservation of angular 
momentum in the particle creation process, the expectation value of the kinetic 
helicity must vanish. However the r.m.s. value, or variance, can be different from 
zero, and this is what we show in this subsection. 

From equation (|3.20| ) we can write a vorticity operator as 



w 



{Ok} 



and define a kinetic helicity operator as 



(P + P) 



{0i}o rp{0k} 



(3.37) 



which in terms of the fields reads 
H 9 



, rp\Wi p. rp 

$ ^ijk 1 - <J> u 3 ± i 



(3.38) 



m 



(p + pY 



ijk 



+ $ (di&) (9,$) (dfl) + (<9 fc $) (a,-$t) 6 (d&) 



(3.39) 



Observe that in principle it does not vanish identically because there are two fields 
involved in its expression. The r.m.s. value of Hf is again given by the vacuum expec- 
tation value of the Hadamard two point function, i.e., (0 | {Hf (f, r) , Hf (f', r')} | 0) 
from where we obtain a spatial two point function as 



Ef (r, x, x') 



lim (0\{Hf(x,r),Hf(x',r')}\0) 



(3.40) 



The calculations are rather long but straightforward; details are given in Appendix 
D. We quote here the main results. When replacing the fields we obtain, as in the 
case of the velocity correlation Rij, several kernels: one with vacuum contributions 
only, another with mixed contributions from vacuum and created particles, and a 
third with contributions from the created particles only. Terms containing \(3 k \ 2 give 
the main contribution, because, as was the case for R^, terms with a K a* m . . ., etc. 
oscillate, and will give negligible contributions when integrated. In Fig. (|I|) we show 
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Figure 1: 16-7T 6 [a (r) /H] 2, Ef (t, £) as a function of £, for fixed m/H = 10 6 and tq - 
0.7 x 10~ 9 (large dashing), tq = 10~ 9 (medium dashing) and tq = 1.3 x 10~ 9 (tiny dashing 
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Figure 2: 167T 6 [a (r) /H] 2 Zf (t, |) as a function of £, for fixed to = 10~ 9 and m/il = 
0.8 x 10~ 6 (large dashing), m/H = 10~ 6 (medium dashing) and m/H = 1.2 x 10~ 9 (tiny 
dashing) 

the dependence of 167r 6 [a (r) /if] 2 5* (r, £) on £ for fixed m/H = 10~ 6 and three 
values of r and in Fig. (H) the dependence on £ for fixed r = 1CT 9 and three values 
of m/H. 

We can estimate the (dimensionless) coherence length of the kinetic helicity as 



(3.41) 
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In Fig. (|^) we show A 2 as a function of m/H for To fixed, and in Fig. and (f|) the 
converse. From both figures we see that, for the chosen values of the parameters, 
A ~ 1, i.e. it is of the order of the particle horizon's scale. 
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Figure 3: A 2 as a function of tq, and for m/H = 0.8 x 10 6 (large dashing), m/H = 10 
(medium dashing) and m/H = 1.2 x 10 -9 (tiny dashing) 
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Figure 4: A 2 as a function of m/H, and for tq = 0.7 x 10 9 (large dashing), tq = 10 9 
(medium dashing) and To = 1.3 x 10~ 9 (tiny dashing) 

The r.m.s. value of the kinetic helicity is obtained by volume averaging expression 
( p.40|) over x and x', i.e. 

S L = [ d [vol (X)] — i- f d [vol (x 1 )} H* (r, x, x') (3.42) 

V Ol \X) Jvol(x) V OL{X ) Jvol(x') 

Kinetic helicity is a global quantity, that can depend at most on the dimensionless 
characteristic scale C < 1 of the integration volume. To evaluate the integrals in eq. 
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( |3.42|) we can proceed as follows. As we are considering scales < 1 we can develop 
eq. ( |3.40|) in Taylor series to second order in x — x' 1 , then using eq. ( 3.25|) we have 



Z Hc (r,£) 



8Hrl 
7ra (t) 



3 /iP 2 
2 \m 



9 



1 



16/ 24 
C 2 



240 A 2 (m/H, r ) 



1/2 



where C = HL, and with 



(3.43) 



2187 1 
" 2097152 



x 



x 



1 H 



H 

m 

2187 1 



1 (E 



2 \ m 



+ 



9477 



29360128 r 2 



1 2 



12 m 5242880 r 2 



3 fH 

8 \ m 



+ 2\m 



+ 



3 Z-ffV 2187 1 
8 \m) + 4194304^f 

9477 1 1 ^ 



29360128 r 2 



the zeroth order in the Taylor expansion of 5* (r, Q . 



(3.44) 



4. Characterizing the turbulence: viscosity, Reynolds number 
and characteristic decay and correlation times 

Turbulence sets in at the inflation-reheating transition, and afterward, as there is 
no stirring forces acting on the flow, it decays. Since the fluid is made of particle - 
antiparticle pairs, annihilation also occurs. We have two competing processes, whose 
characteristic times must be compared in order to decide which one dominates: the 
decay time of turbulence and the time of particle anti-particle annihilation. 

To properly characterize the turbulence, we must first calculate the Reynolds 
number of the flow. This number is defined as Re = ul/v, with u and / characteristics 
velocity and scale respectively and v the kinematic viscosity. As it may happen that 
turbulence is not fully developed at all scales, we must calculate this number for each 
scale. The scale dependent Reynolds number can be written as 

i2e(«) = ^I^l (4.1) 
v 

where v = r\j (p + p) is the dimensionless ratio of the fluid shear viscosity to the 
energy density k^ 1 = a (r) H^ 1 ^ 1 is the scale of interest and v K ~ a/ 2E k (r) the 
estimate of the velocity at the corresponding scale. To estimate rj we follow the work 
of Son and collaborators p5| , and consider it proportional to the entropy density, 

1 It can be seen from eq. (p.6|) that this is the next to leading order. 
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i.e., rj/s — 1/47T. We take s as proportional to the quasiparticle number density n, 
i.e., s ~ n (r), with n (r) 



# 3 1 
n (r) ~ — 



Q> [ T ) Tn 



3H 3 2 /9\ 4 vr 



2m 3T ° + 1 16 J 28 



(4.2) 



(see Appendix B). Using eq. ( p.25[ ) we have that 

_ 1 a (r) [(3/2) (H/mf r 2 + (9/16) 4 vr/28] 
U (2vr) 1/2 # [(3/2) (H/mf r» + (9/16) 4 /24] 1 ' J 

Replacing everything into eq. (f4.1|) we obtain 

Re U r) * (2,)^ ( ^] " ,/X ^ ^ (gM Z° JL (4 4) 

1 16^ V 256 [(3/2) ( J ff/m) 3 r 2 + (9/16) 4 vr/28] 

We see that, for k ~ 1, we can have Re ^> 1 if 



i.e., if the transition between inflation and reheating is very fast. In Fig. (H) we 
plot Re (k, t) /a 1 / 2 (r) as a function of n for fixed m/H and three values of r , it is 
seen that Re increases as the duration of the transition inflation/reheating shortens. 
In Fig. (^) we plot the same as in the first figure, but with fixed r = 10~ 9 and 
m/H = 10 -5 , 10~ 6 and 10 -7 , and we observe that Re diminishes with decreasing 
m/H. In both figures Re is a decreasing function of k, hence only the modes near the 
horizon can be considered as turbulent. In Fig. (0) we plot Re (k = 1, r) /a 1 / 2 (r) as 
a function of r and for m/H = 10~ 5 , 10 -6 and 10~ 7 . In this case Re has a peak at 
certain value of r , and this peak is higher and occurs at shorter values of r as m/H 
diminishes. As tq grows Re decreases monotonically. Finally, in Fig. (|8p we show 
Re/a 1 / 2 (r) as a function of m/H, for k — 1 and tq = 10~ 7 , 10~ 8 e 10 -9 . Again Re 
peaks at certain values of m/H and the peak is higher and occurs at smaller values 
of m/H as r diminishes. 

The decay time of each turbulent mode is given by td (k) = 1/vk 2 = a 2 (r) H~ 2 /vk 2 . 
Using eq. ( |4.3|) we can write 

(2nf 2 a(r) [(3/2) (H/mf rj + (9/16) 4 /24] 1 
d [K) ~ H [(3/2) (H/mf r 2 + (9/16) 4 tt/28] k 2 1 ' 

To estimate the pair annihilation time of each mode, we consider that the par- 
ticles are relativistic, according to the result obtained for their energy density. On a 
dimensional basis, this time can be estimated as 

ta(K)~—^ (4.7) 

n (t) a K u r 
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Figure 5: Re (k, r) /a 1 / 2 (r) as a function of «, for fixed m/H = 10 5 and to = 10 7 (large 
dashing), 10~ 8 (medium dashing) and 10 -9 (tiny dashing). Re increases as tq diminishes 
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Figure 6: Re(K,r) /a 1 / 2 (r) as a function of k, for fixed tq = 10~ 9 and m/H = 10 -5 
(large dashing), 10 -6 (medium dashing) and 10~ 7 (tiny dashing). Re diminishes as m/H 
diminishes 



with n (r) the particle density, a K the annihilation cross-section and u r the relative 
velocity between species which we take u r ~ 1. We make a crude estimation of 
the cross section, as being the same as for e + e~ annihilation for 7 ^> 1. i.e., 
a — 7r7"o/7, with, r = a/m, a = 1/137 the fine-structure constant, 7 = e/m with m 
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Figure 7: Re(K = l,r) /a 1/2 (r) MS M function of To, and for m/H = 10 5 (large dashing), 
1CP 6 (medium dashing) and 10 (tiny dashing). Observe that Re peaks at a certain value 
of To- The peak is higher and occurs at smaller values of tq as m/H diminishes. 
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t igure 8: Re(n = l,r) /a 1 / 2 (r) SIS Si function of m/H, and for tq = 10 7 (large dashing), 
10~ 8 (medium dashing) and 10 -9 (tiny dashing). Observe that Re again peaks at a certain 
value of m/H. The peak is higher and occurs at smaller values of m/H as tq diminishes. 



the particles mass and e ~ HVt K = H \k 2 /a 2 (r) + (m/H) 2 ^ ^ the maximum energy 
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of each mode. We then have 



mH 

replacing in eq. (|4.7|) we have 



a 2 a (t) 

°« - 1 = ( 4 - 

k 2 + a 2 (r) (m/H) 2 



a 2 (t) mTn 



ta (*) ~ -dr-sf r^!., ^3 2 , ,,w^4 (4-9) 



^//t 2 + a 2 (r) (m/H) 
vra 2 if 2 [(3/2)(F/m) 3 r 2 + (9/16) 4 /28] 

Comparing both times we have 

t d (n) (2vr 3 ) 1/2 a 2 H [(3/2) (H/m) 2 r 2 + (9/16) 4 /24] 



(«) a(r) mr K 2 y J K 2 + a 2 ( r ) (m/tf) 5 



(4.10) 



as if/m 3> 1 and r 2 C 1 we have that td (k) /t a (n) 1, i.e., annihilation occurs 
before the end of turbulence. 

It can be seen that t a is also much shorter than other time scales pertaining to 
the flow, such as the ratio between the radius of the largest turbulent eddy (i.e., the 
horizon as it is there where Re (k) takes its largest value) to the velocity associated 
to that scale. Therefore in what follows we take t corr = t a . 



5. Magnetic field amplification due to dynamo action 



According to eq. ( |2.8| ) we must now evaluate the amplification exponent, Tfhi 2 (r, C) t (k) /La (r), 
where La (r) is the physical coherence scale of the magnetic field and £ is the (dimen- 
sionless) coherent scale of the kinetic helicity. From the discussion in Section III, the 
shortest time is the annihilation time, which for scales such that k ~ 1/C > m/H 
reads 



t a (£) 



a T TflTft 



1 



na 2 H 2 [(3/2)(#/m) 3 r 2 + (9/16) 4 /28] £ 
Writing L = if _1 £ the amplification exponent reads 



(5.1) 



£ Wc (r,£)t a (£) 8 



La It) 



7r 3 a 4 



a 3 (r) 



2^8 



m t, 



o 



H 2 [(3/2) (H/mf t 2 + (9/16) 4 1/24] : 
1 



x ■ 



[(3/2) (H/mf r$ + (9 /16) 4 /28]' 



1 - 



c 2 



nl/2 



240 A 2 (m/H, r ) 
which by simple inspection is seen to be very small. 



(5.2) 



- 18 - 



6. Discussion and Conclusions 



In this paper we have studied the possibility of turbulent dynamo action during 
reheating. We considered the presence of a charged scalar field minimally coupled 
to gravity. This field is in its invariant vacuum state during inflation. When the 
transition to reheating takes place the vacuum state turns into a many-particle state. 
For sub-horizon modes of the field, the number of created modes depends on the 
details of the transition. Therefore during reheating, besides the decay products 
of the inflaton we also have a plasma of scalar particles which is at rest in the 
comoving frame. We characterize the fluctuating velocities of this plasma giving 
their spatial two point correlation function and the kinetic energy associated to each 
Fourier mode of the stochastic velocity field, eq. (|3.34|) . We evaluate the Reynolds 



number associated to each mode, Re (k), which turns out to depend on the physical 
parameters of the problem, namely m, H and the duration tq of the transition 
from inflation to reheating. If r is small enough, then there is a range of k for 
which Re (k) > 1 and the flow can be considered as (mildly) turbulent. As there 
are no stirring forces, the turbulence we refer to decays, each mode doing so in a 
characteristic time td (k) given by eq. ( |4.6| ). Besides as the plasma is a particle 
anti-particle one, each mode of the scalar field (not to be confused with modes of 
the stochastic velocity field) annihilates in a characteristic time t a (k) given by eq. 
( [4.9|) . When comparing both times we find that annihilation dominates over decay, 
eq. ( f4.10| ) and hence for practical purposes we can consider the turbulence as steady. 

The sufficient condition to have a large scale kinematic dynamo is the flow to 
be endowed with kinetic helicity |35 |. The non-trivial result of this paper is that 



the scalar plasma does possess a non null rms kinetic helicityeq. ( |3.43|) From Figs. 



([!]) and (|j) we see that, for the parameters for which Re (k) > 1, the characteristic 
scale of the kinetic helicity is of the order of the particle horizon, thus allowing for 
kinematic dynamo action. 

The existence of an rms helicity is due to the presence of the two scalar fields, 
$ and as is evident from eq. ( |3.39| ). Moreover, though the helicity may have 
either sign, in the average the amplification effect dominates. From the simplest 
model of kinematic dynamo, eq. (|2.7|) , we compute the amplification factor of an 
initial seed field, eq. ( j5.2|) , and find that for the physical parameters of the scenario 
considered, it is very small. In spite of this result, we believe our work shows the 
need for exploring the impact of nonlinear effects in the early universe. These effects 
offer the most natural answer to the riddle of the survival of the primordial magnetic 
field until the epoch of structure formation, in spite of the 1/a 2 damping induced by 
the Hubble expansion. 
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A. Bogoliubov coefficients 

We assume that during the reheating period the scale factor of the Universe scales 
as t 2//3 |^9|, |5(|, while for inflation we consider a spatially flat de Sitter universe. 
For large wavenumbers the Bogoliubov coefficients are sensitive to the details of the 
transition, while for small wavenumbers the coefficients can be found assuming an 
instantaneous transition. This dependence on the transition details for subhorizon 
modes was also recently analyzed by Zaballa and Sasaki |HJ] in the context of creation 
of metric perturbations at the end of inflation. 

The Klein Gordon equation for a free field in a FRW Universe is 



k 2 /m\ 2 3 



a 



a a 
2 \ a + 2a 2 



It is seen that 



77 - + 



1 rfV/ 2 _ 9 
2a 2 ) ~ a 3 / 2 dr 2 ~ 4* ^ 



(a.i; 



(A.2) 



with x = (r — T\) /t where r is the time the transition lasts. We assume Tq T\ -C 1. 
In terms of x eq. ( |A.2|) gives 



d 2 o?l 2 9 
dx 2 



^a^fix) 



which can be integrated giving 



a 3 / 2 (x) 



A + Bx + -r 2 
4 u 



dy (x-y)f (y) a 3/2 (y) 



(A.3) 



(A.4) 



where x\ = Ti/tq. The constants of integration are obtained by matching to the 
inflationary solution at r = 0. We get A = 1 + tqXi and B = tq, then 



a 3 / 2 (x) 



with 



9 f x 

1+to {x x + x) + -r 2 / dy{x-y)f (y) a 3/2 (y) = F (x)+r xG i 

^ J —X\ 

F (x) = 1 + r oXl - jT 2 [ dyyf (y) a 3 ' 2 (y) 

^ J —XI 



G(x) 



, 9 



dyf (y) a 3 / 2 (y) 



x) (A.5) 

(A.6) 
(A.7) 
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For x > 1 we have / ~ 0, hence F and G are constants in that x range. We now 
write K.G. equation (|A.1|) as 

d 2 ^ K 



dx 2 



+ r 2 ^ (x) ij K = 



with 



m 



9 



(A.8) 



(A.9) 



We are interested in the behavior of the solutions equivalent to Hy (x) (the positive 
frequency solutions for a spatially flat de Sitter spacetime) in x = —x\. There are 
two possible situations. Given that x = is the middle of the transition, we consider: 
(a) Q (0) t < 1: the details of the transition are not important; (b) Q (0) t > 1: 
the details of the transition matter. For modes inside the horizon (k > 1) we can 
consider the WKB solution 

giro S [a:] 



V«+ (z) 



with fl (x) = dS [x] /dx. The derivatives are 



(A.10) 



d 2 ip K+ 
dx 2 



d^ K ^ 
dx 



ir {l K (x) + 



1 K (x) 

2 Vl K (x) 



KM 



1 fWJx) 



4 \njx) 



l fn'Ax] 

2 



and then the equation for -?/> K reads 



dx 2 



(X) 



i /fiUaO 



4 \njx) 



K fr) 

Q K (x) 

2 1 

+ 2 



(A.ll) 
(A.12) 



2 

^(x 



£l K (x) 



(A.13) 



The solution can be expressed as a superposition of positive and negative frequency 
modes as 



ip K+ (x) H (p K+ (x) 

TO 



XI 



dyip K - (y) 



i fn' 



4 \n K 



+ 



TO 



<Pk-(x) / dy^K+iy) 



+ 



i (si' 



2 \n K 



f2 K 



^ic(2/) 



(A.14) 



Whereby we read the Bogoliubov coefficients 

z r i fo! 

a K =l + — / dy<p K -(y) — -± 
T~0 J - Xl 4 V i2 « 



i /n: 



fa (y) 



(A.15) 
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T 



Pk = / dyip K+ (y) 



X\ 



i fn'\ 2 i /n 



4 \n 



2 V tt 



tPAv) (A.16) 



which should satisfy |a K | 2 — |/3 K | 2 = 1. To obtain a simpler expression, we consider 
an iterative solution. To lowest order, i.e., ip (y) ~ <p + (y), we have 



a 



(°) = 1 + 



T 



J-i 



dyp K - (y) (f K + (y) 



■n 



i fmv i 

+ 



4\fl 



2 I ft 



/3f - -- r dyvU (y) 



4 UJ + 2 U« 



and 



^ K (x) = 4 V+ (a:) + 4°V«- (*) 
Performing another iteration we obtain 



a^. 1 ' ~ exp 



r 



— / dyip K+ (y)(p K -(y)Q K (y) 



■t-i 



1 + 



(0)1 2 



(A.17) 
(A.18) 
(A.19) 

(A.20) 



Observe that it is not necessary to perform another iteration for (3^. The normaliza- 

2 2 (0) 

tion condition \a K \ — \(3 K \ = 1 is satisfied up to a term (3 K /4, indicating that 



this coefficient must be 



(0) 



<C 1 in order to render our expressions valid. 
Integrating by parts in eq eq. (|A.16|) and neglecting surface terms we have 

1 A£Y 



dy-^v*+ (y) <p' k + (v) 



(A.2i; 



-ZTnfL 

2 0, 



Observe that the first term is suppressed by a factor of r , so we shall not consider 
it further. We are now at the point where details begin to matter. Write 

-2k 2 a' 9,. " 

-J ( x ) 



a 2 a 4" 



(A.22) 



The k 2 term may be neglected even when k 2 is large. To see this, observe that 



HJ 4* 



so 



n K 



2 U /(l) -(«) )-a + l f(X \ 



Writing 



2^ 

2 (a 3 / 2 )' _ 2 t G 
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a 3 a 3 / 2 3 (F (x) + t xG {x)) 



(A.23) 



(A.24) 



(A.25) 
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we see that this function is suppressed by r . Therefore in eq. (|A.24 ) the only term 
that is not suppressed is the last one. So we finally have 



ft? 



(°) 



T 



9 

16 



dy 



f(y) 
IW (y) 



(y) 



(A.26) 



We see that /3 K is essentially the Fourier transform of (/') . Since this function 
has a peak whose width is ~ 1, by Heisenberg's principle we expect to get a negligible 
result for r Q K ^> 1, namely for k ^> r^ 1 . Observe however that this scale can be 
extremely high. To give concrete results, let us consider /' = const = —1/2 and 
assume that we can make a linear approximation in the exponent of eq. ( |A.10| ), 
S [x] ~ S [0] + Q [0] x. Assuming x\ ~ 1 and that fl K (y) is a slowly varying function 
of time to keep only the surface terms in the integral, we obtain eq. ( |3.16| ). We see 
that the number of created particles with large k is very sensitive to the details of 
the transition between inflation and reheating; it would actually diverge in the limit 
ro i — y 0, which is therefore unphysical. 

For small k an instantaneous transition can be considered, and the coefficients 
calculated by directly matching the inflationary and reheating solutions at r = 0. 
Assuming again a WKB form for the modes during reheating and the usual Hankel 
function for de Sitter 



. For this transition the full expression for j3 K is 



7T 



1/2 



vmH 2 (o) 



kH, 
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v-l 



+ HP 



L3n* (o) 



— v — iQ K (0) 



(A.27) 



with v = ^9/4 -m 2 /H 2 ~ 3/2 - m 2 /3H 2 . For k < 1 the Hankel functions can be 
approximated as 



Using 



DT(u-l) 



7T ' ' \2, 

and the fact that 2k 2 /VL 2 k (0) 3 - v 



( K ) ~ --Y{v 



(A.28) 
1 we get eq. 



B. Calculation of (p + p) and n 

In terms of the scalar field we have 

, , v „ 4 gg) + (Tj a ) 
{p + P) ~ ^ 

Using eq. ( p.3|) we obtain 



3 ) + — (V$ • V$ f > 



Replacing the decompositions ( |3.10|) and using 

[a K ,al] =(2nf 2 a 3 (T)5(R-zu) 



(b.i; 



(B.2) 



(B.3) 
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we obtain 



— <V$- V$ f ) 
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# (r) (r) + 0i (r) (r 
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(B.4) 
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a 2 (r) 



Using decomposition ( |3. 15|) we identify two different contributions to the integrand: 
one from pure vacuum 
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and one from the created particles, 
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We are interested in M% and of it, the contribution from the \f3kf terms is the most 
important: as a k oscillates (see eq, | |A.20| |), the terms proportional to a k and a\ will 
give negligible contributions when integrated. Replacing the WKB form for (r), 
eq. ( p.!4|) we obtain 
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Replacing 



a a 



3/2 : 



Oc (r) 
ttk (r) 



1 7 / 2 a. 



"a^n»(r) - 



a 7 ^ (r) 



(B.8) 



we have that 
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By simple inspection we can see that the terms that contribute the most are 
3 /$$t\ + J_(V$- V$ f ) ~ 
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because they decay more slowly than the others. We must now replace the Bogoliubov 
coefficients and perform the integrations. For long wavelengths, i.e., those in the the 
interval (0, 1) we use the expression (|3.18|) . Thus in this case we must evaluate 



(o) 



(^-i) 2 r 2 (i/2) 

(2vr) 3/2 7ra 3 (r) 



dn 



V K (0) K 2v 



Sn k (r) + 



a 2 (t) n k (t) 



(B.ll) 



Since m/ff < 1 and we are considering a period of time in which a (r) does not 
differ very much from unity, we can take a (r) ~ 1 in all the roots, and so we have 
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And finally the contribution from long wavelengths reads 
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To evaluate the contribution from the short wavelengths we use eq. ( 3.17 ) in eq. 
( [B.10D , so we have to compute 
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As in this case m/H <C k, we can neglect that term and so the contribution from 
short wavelengths reads 
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Gathering expressions ( |B.13| ) and ( |B.15| ), approximating u — 1 ~ 1/2 and using that 
r (1/2) = 7T 1 / 2 we arrive at expression fl3.25| ). 

For n, the number density of created particles, we have [EM |39|, HO 
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Using again eqs. ( 3.18 ) and (ft. 17 ) in the appropriate momentum intervals we obtain 
eq. (Ob. 



C. Calculation of the velocity correlation spectrum 



We start by replacing eq. (3.23) into (3.22) and noting that three kernels build the 
correlation function: one with the vacuum contributions 
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another with mixed contributions from the vacuum and the created particles, 
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and a third kernel with contributions from the created particles, 
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where the dots in square brackets indicate more terms with combinations of a K , (3*, 
a-,;,, filj. Of the three kernels, ( K ' w i r ' r ') gi yes the main contribution, because 
it has no vacuum contribution. Observe that as the coefficient a K is oscillatory (see 
eq. | A.20 ] in Appendix A), the terms with coefficients with a K and a m will give 
negligible contributions when integrated. Therefore in what follows we shall analyze 
only the terms in |/3 K | 2 |/3 ro | 2 . From direct inspection of eq. ( |C.3|) we see that the 
only terms that can survive after integrating are those proportional to vj l vji and to 
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k 1 k? . Thus we have to evaluate 
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The modes during reheating are of the WKB form, and thus the G^ + (t, r') reads 
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The velocity spectrum is defined in eq. ( |3.26|) , so taking the coincidence limit r = r' 
of eq. ( |3.24| ) and transforming Fourier we obtain 
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where we used the isotropy of the Bogoliubov coefficients to replace 

«V ^ ^ 2 5 lj (C.7) 

as those terms are the ones that give non null contributions. After replacing the 
propagators and their derivatives, the velocity correlation can be written as 
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Here we must replace the Bogoliubov coefficients eqs. f p.171) and ( |3.18| ). We are 
interested in short wavelength velocity modes, i.e., those inside the particle horizon 
for which q > 1. However, care must be taken when ac approaches q as in this case the 
Bogoliubov for long wavelengths must be used. After long but straightforward cal- 
culations we obtain that the full expressions for the contribution of long wavelengths 
to the velocity spectrum is 
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while for short wavelengths we have 
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As ^ > 1 the main contribution comes from term whose inverse power of ? is the 
smallest, so we shall keep only them. Observe that they come from the contribution 
of short wavelengths, and therefore depends strongly on the details of the transition 
inflation-reheating. The two contributions to the velocity spectrum are then 
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Replacing eq. ( |3.25|) we obtain eq. ( |3.27|) . 
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D. Calculation of the kinetic helicity 



We start from eq. ( ft.40|) and replace decomposition (|3.10|) of the fields, thus obtaining 
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Noting that to avoid an odd integrand, we must contract a' with a, we are left with 
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where we used [a K , Considering all possible combina- 

tions of the remaining moments, we finally have 
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with £ = x — x' . Replacing 0£ (r) = a fe 0f (r) + (5k4>k * ( r ) we obtain, as in the case of 
the velocity correlation R^, several kernels: one with only the vacuum contribution, 
another with mixed contributions from vacuum and from the created particles, and 
a third one with the contribution of only the created particles. The expressions are 
rather long, but they are straightforwardly obtained . Of the one due to the created 
particles, the part with \(3k\ 2 gives the main contribution, because as was the case for 
Rij, terms with a K a^ . . ., etc. oscillate, and will give negligible contributions when 
integrated. Therefore we shall consider 
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Replacing the WKB form for the modes and keeping only the slowly decaying terms 
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we can express eq. ( D.4|) as 
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Working in spherical coordinates and performing the angular integrals for each mode 
we are left with 
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These integrals can be performed straighforwardly using the same approximations 
as for obtaining 
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and 
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